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SENSOR FOR MEASURING INSTANTANEOUS ANGLE OF ATTACK 
OF HELICOPTER BLADES 


By 

P. Stephen Barna 1 and Henry W. Liu 2 


SUMMARY 


During the period July 1976 to July 1977, the various research activities 
may be divided into three different categories: 

1. improving and testing probes; 

2. theoretical studies of probe motion; 

3. improving research facilities. 

Substantial progress was made in all aspects and it is anticipated that a 
satisfactory solution to the problem of measuring angle of attack of helicopter 
blades may be found in the near future. 

SYMBOLS 


b 

probe constant for a real fluid flow 



b o 

probe constant for ideal flow, b Q = 2.25 



bi 2 , b23 

probe constants based on portholes 1, 2 and 

2, 5, respectively 

b 

r 

the ratio of b 2 3 and bi 2 , b^ = b 23 /b 12 



e l 2 j e 23 

functions that represent the dynamic effect 
AP2 3 

for Ap 12 

and 

f 

frequency, Hs 
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Subscripts 

i 

o 

U 

L 

D 


gravitational acceleration, 9.81 m/sec 2 

static pressures at portholes 1, 2, and 3, respectively 

static pressure at point i on a sphere 

total pressure 

free stream static pressure 

pressure differences between portholes 1, 2 and 2,5, respectively 
period of one complete oscillation (T = 1/f) 
free stream velocity, m/sec 
time, seconds 

velocity components in x, y, and z directions respectively 

air velocity at point i on a sphere 

probe velocity under oscillating motion, m/sec 

relative velocity at probe upstroke and downstroke respectively 

stagnation point 

density 

pitch angle, degrees 
viscosity 

reference displacement angle, degrees 
half of the amplitude angle, degrees 
pitch angle at upper and lower probe positions 
velocity potential 

input 

output 

upper 

lower 

displacement 


u 


upstroke direction 
downstroke direction 


d 

t total 


DESCRIPTION OF IMPROVED PROBE 


The first, now outdated probe, referred to as "ODU #1,” was a three-porthole 
spherical probe where the two transducers were arranged symmetrically as shown 
schematically in figure la. In this arrangement the pressure differences between 
ports 1 and 2 and 1 and 3 were measured across the transducers T^ and T 2 
respectively, where port 1 was common to both. Results of experiments performed 
on probe 1 were reported in Technical Report 76-T14, July 1976. 

The second probe referred to as M ODU #2" is also a three-porthole spherical 
probe of the same dimensions as SI, where the pressure difference between ports 

1 and 2 was measured with transducer Tj , while the difference between ports 

2 and 5 was measured with transducer T 2 , is shown schematically in figure lb. 
Details of the probe are shown in figure 1c. 


According to the theory of flow around a sphere, direct measurements of 
pressure difference p 2 - P 3 and p^ - p 2 lead to the angle of attack, e. In 
the July 1976 report, it was shown that 0 can be calculated with equation C7) 
of that report as 


8 = j tan -1 


2 CPl - P2) 


(p2 ~ P3) 


Subsequent studies indicate, however, that the probe constant b Q , calcu- 
lated on the basis of ideal fluid flow around the sphere, is somewhat different 
for real fluids. More particularly, the constant b 2 3 for the pressure dif- 
ference p 2 - P 3 was found to differ from the constant b ^ 2 related to the 
pressure difference p^ - p 2 . Theoretical considerations presented in this 
report show that the angle of attack can be more accurately calculated with 
equation & of Appendix Ai 


f 


P2 " P3 

CP2 - P3) + 2b r (p ! - p 2 j 

where the ratio fa^ = b 2 3 /bi 2 is to be established from calibration experiments 
performed under steady flow conditions. 

Three probes were tested during the period July 1976 to July 1977: ODU #1, 

ODU # 2, and a probe designed by the Flight Instrumentation Division of NASA/ 
Langley Research Center. As indicated earlier, #2 superseded #1, while the 
first NASA probe proved unsatisfactory on account of the misalignment of the 
central port. This fault was corrected and the improved probe proved satis- 
factory. Only the results of probe ODU “2 are presented in this report because 
of space limitations, 

IMPROVEMENTS TO THE TESTING EQUIPMENT 

Improvements were made in two directions : (a} to produce a truly sinusoidal 

oscillating motion, and (b) to increase the air speed in the tunnel. 

a. The first mechanism producing oscillatory motion was described in 
Appendix A of the July 1976 report. This mechanism was subsequently replaced 
in 1977 by a second mechanism called "Scotch Yoke," which is known to produce 
sinusoidal motion. In a "Scotch Yoke" the circular motion of a pivot is trans- 
formed into a reciprocating motion. The pivot is inserted into an elongated 
slot cut into a crosshead that is provided with guides. As the pivot rotates 
the slot translates parallel with itself. Details of the Scotch Yoke are 
shown in figure 2a. The mechanism was designed to fit the existing experimental 
apparatus. 

b. Improvement to the wind tunnel was necessary because of the inefficient 
fan and unsatisfactory "Vee belt" drive. Access to the pulley to the drive shaft 
was attained through a large, 14-in. diameter, funnel, which partially obstructed 
the flow downstream from the fan blades. Subsequently, the fan blades were 
redesigned and the funnel was eliminated by changing the "Vee belt" drive to 

a front-end direct-drive through a floating shaft. The new fan was provided 
with a redesigned set of downstream guide vanes which were housed inside a 
new and much improved housing [fig. 2b). The process of redesigning, together 
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with the problems encountered in the process of manufacturing the components, 
plus the installation, took longer than a year and entailed a massive coordina- 
tion of effort to complete. ' The performance of a new fan proved 100 percent 
satisfactory as it increased the air speed from about 65 ft/sec to 100 ft/sec, 
and the dynamic pressure, q, by a factor of two 1 

THEORETICAL STUDIES OF PROBE MOTION 


The probe motion at present is still subject to theoretical studies. The 
oscillating motion of the sphere on a semicircular path may be considered as 
the result of three component motions: up and down (vertical component); 

back and forth (horizontal component); and turning about its horizontal axis 
(rotary component) . Solution to the horizontal motion has been found in the 
literature, and a similar solution to the vertical component may be found in 
the near future. The rotary component will probably cause no great problem. 
However, the difference between theoretical (potential flow) solutions may 
differ from the real effects, and comparative experiments may become necessary 
on the component motions to verify theory, 

EXPERIMENTAL RESULTS 
Steady Flow Tests 

Results of steady-flow tests are based on the pitch angle 8 measured from 
the horizontal datum line. The stagnation point S shown in figure 5a is only 
applicable to steady flow and is located where the horizontal velocity vector 
U contacts the surface of the sphere. In the steady-flow test, the pitch 
angle was varied from the lowest, 6^ = -30, to the highest (upper), 0y = +50, 
by steps with the probe "fixed" in space while the observations were made. 

Tests on the ODU #2 probe were performed under various air speeds in the 
range of Reynolds numbers* 2.37 to 6.25 x 10 4 (corresponding to air velocities 
of 58 to 100 ft/sec) . In figure 4 the normalized pressure difference Ap 2 3 /q 
and Api 2 /q are plotted against pitch angle 8, More particularly in figure 


* The Reynolds number may be calculated from the formula Re = 625 U, where U 
is in ft/sec. 
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4a both pressures are plotted against 0, while in figures 4b and 4c the same 
quantities are plotted against sin 20 and 1//2 sin (tt/ 4 - 20), respectively. It 
appears from figure 4b that the probe constant b 2 3 = 1.8 to 1.9 and from figure 4c 
b 12 = 1*9 to 2.0, the exact value depending on the Reynolds number used. The 
ratio b^ = ^2l^\2 is plotted against the Reynolds number in figure 4d. 

From the results thus obtained it is possible to predict the output angle 
using equation (8) given in Appendix A. In figure 5 the output angle 
6 out is plotted against the input angle 6^ n , the latter being measured 
with a clinometer, while the former was calculated assuming b r = 0.857. 

Some disagreement with 0^ is noticeable due to experimental error. The 
discrepancy, however, may be considered small. Variation in b 2 3 and* b^ 2 
were found to be strongly dependent on Reynolds number, especially in the 
subcritical range (Re < 2 x 10 5 ) ; no data in the supercritical range 
(Re > 2 x 10 5 ) are available at this time. 

Tests With Oscillating Probe 

Tests with the probe oscillating fall into the unsteady flow category 
because the flow pattern about the sphere changes with time. Under oscillating 
motion the spherical probe head follows a circular arc path enclosing an ampli- 
tude angle 2a. Measured from the datum are the central "extreme" angles 
and 8 U , at which the sphere velocity, -v t , is zero; while at the displace- 
ment angle 6^, the sphere moves with maximum speed. It may be seen from 
figure 5a that 


1 a 


0 L + S U 


and 


0_ = a - 0. = 0 


U 


- a 


The probes were tested in a frequency range 1 to 5 with air speed varying 
from 58 to 100 ft/sec. Amplitude was set at a nominal 2a = 18 degrees, while 
displacement, 0 D> was varied from 0 to 15 degrees. 

Results of tests performed under oscillations are presented in figures 6 
to 9, where the normalised pressure differences Ap 23 /q } Apj^Al are plotted 


against 0. The abscissa consists of two parts: the upstroke and the downstroke 

The upstroke motion starts with 9 = 0^, shown at far left, and finishes at 
0 = 9 jj at center; the downstroke starts with 0 = 0y at center and finishes with 

0=9 at the far right. The time for the complete period is T = 1/f. 

Li 

The difference in the curves between the up and down stroke may be partly 
explained by first considering the velocity vector diagram shown in figure 3b. 
During the upstroke, the stream velocity U and probe velocity -v combine 
into the relative velocity V , while during the downstroke they combine 
into V^. When displacement 0^ is different from zero and 

V u >:> V d’ hence t * ie Pressure differences change accordingly. In addition to 
the change in relative velocity the stagnation point also changes location. 

During the upstroke S' will be located below S (S being the steady flow 
stagnation), while during the downstroke S" will be located above S.* 

Further, theory indicates that under unsteady flow conditions the steady 
flow Bernoulli equation must be corrected by adding the term d^/dt, where $ 
is the velocity potential. Accordingly, the pressure difference between ports 
increases above the value calculated for steady flow conditions when f = 0. 

The curves presented in figures 6 to 9 are results of "parametric” studies. 
Because of the variables U, f, 2a, 0 Q , and probe geometry, it is only 

possible to study the effects by varying only one parameter at a time while 
keeping all others constant. During one set of experiments, the probe geometry, 
2a, the Reynolds number of the sphere, and 9^ were kept constant leaving 
only f as variable. The experiments were repeated for various values of 0^. 
Still, in another set of experiments f was kept constant and the Reynolds 
number was varied, and so on. The set of data thus obtained were plotted on 
each diagram for a particular Re and 0 , where 2a = 18 degrees for all 
figures presented (figs. 6, 7, 8, and 9). Note that if 2a were to be 
replaced by the time , T, of a complete cycle of oscillations, the abscissa 
spreading over 2a would inversely decrease with increasing frequency 
(T = l/f) . These experiments were repeated for various Reynolds numbers. 

Since the velocity of the oscillating probe is zero at each end of the 
"stroke" (9^, 9y) , the pressure difference should approximately equal the 
steady flow values. It appears, however, that the steady flow values are 
somewhat lower, which may be attributed to inertia effects in the flow. The 
boundary layer moving in the flow direction during an upstroke probably keeps 

*""Note: §"• and S” are not. shown. _ 
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on moving in that direction for a short period of time even after the motion of 
the sphere is stopped, and may retain its motion after the sphere motion is 
reversed from upstroke to downstroke. 

It appears from figures 6 to 9 that for the lower Reynolds numbers, the 
effect of frequency on the pressure difference is greater than at higher 
Reynolds numbers. This may be explained by considering the relative magnitude 
of the velocity vectors U and v^. At higher Reynolds numbers the length 
of U vector increases [see fig. 3b), and the resultant relative velocity V 
is then largely controlled by U. It may be anticipated that under actual 
operating conditions, when the Mach number attains values of 0.5 at the probe 
location on the helicopter blade, the effect of oscillatory motion on the 
overall performance of the probe will become negligible. This assumption, 
however, needs to be confirmed by experiments on the probe with very high 
velocity air streams which are unavailable at present at Old Dominion 
University. 

EFFECT OF REYNOLDS NUMBER AND FREQUENCY 

Variation of &P 2.3 and Api 2 with frequency may be studied in figure 10. 
The curves shown in figure 10a are of sinusoidal character. The solid line 
corresponding to ’’zero" frequency (very slow motion) is, of course, purely 
sinusoidal and was established from "steady flow" results shown in figure 4, 

The increase of frequency to f - 2 results in a higher peak which occurs 
earlier, and a further increase to f - 5 shows similar effect. Since during 
these experiments the air velocity approaching the probe remained constant, 
one concludes that increase in frequency has two effects: (a) increase in 

pressure difference due to the unsteady flow effects; and (b) the occurrence 
of pressure peaks progressively earlier with increasing frequency. This is 
clearly an out-of-phase effect, and is partly due to the stagnation point 
changing its position relative to the steady flow location shown as S in 
figure 3. The other factor is probably due to inertia effects in the boundary 
layer. These effects are less marked in figure 10b. 

Variation of Ap with air speed (corresponding to Reynolds number changes) 
(see fig. 11) shows for constant frequency a reverse trends as compared with 
frequency changes because the increase in air speed brings the higher air speed 
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curves closer to the steady flow (f = 0) case. This result is of great impor- 
tance because it clearly indicates that performance of the probe under very 
high air speed (M ~ 0.5) conditions may be predicted from steady flow experi- 
ments for a known frequency. 

Experiments with high speed air flow (U ~ 300 ft/sec) have recently been 
started, but no results are available at the time of writing this report. 
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APPENDIX A 


THEORETICAL CONSIDERATIONS 


STEADY FLOW 


For steady uniform flow conditions of an ideal fluid passing over a 
perfectly spherical object, the pressure at a point i on the surface is 




( 1 ) 


where p^ and are the static pressure and the relative surface velocity 

respectively at point i on the sphere, and p is the stagnation pressure. 

The pressure coefficient for a spherical object may be derived from 
potential flow theory and is given by 



1 - b sin 2 8 
o 


C2) 


where 0 is the central angle enclosed between point i and the stagnation 
point S (see fig. Al) . Potential flow theory yields a sphere constant b Q 
equal to exactly 2.25. 

For a real fluid flow, it appears convenient to employ a relation similar 
to equation (2] , 



1 - b sin 2 e 


(3) 


where the sphere constant b is expected to differ from b Q ; however, its 
value needs to be established from experiments. Note that equation (5) usually 
holds only for angles of 8 <_ 60 degrees. 

Pressure differences existing between any two points on the sphere may 
be established with equation (3) . Consider the pressures at the side portholes 
2 and 5 as shown in figure 1. Let the pressure at point 2 be 


p 2 = P t - b q sin 2 02 


and at point 5 


P3 = P t ' b q sin2 e 3 


where 


q.ipU* ■ 

The pressure difference dp 23 = P 2 _ P 3 between side ports 2 and 5 is 
obtained by subtraction: 

Ap 2 3 = P 2 “ P 3 = b PCsin 2 0 3 - sin 2 0 2 ) 

Let the central angle enclosed between ports 2 and 5 be exactly 90 degrees, 
0 2 + 63 = tt/ 2 and let 0 = (83 - 8 2 )/2. Upon substitution one obtains 

dp 2 3 = q b 23 sin 20 

where the value of b 23 depends on the location of ports 2 and 3. 

The pressure difference between center port 1 and side port 2 may be 
expressed as 


&P 12 = Pi - P 2 = -p b 12 sin (j - 20 ) 


where the sphere constant bi 2 depends on the location of the side port 
(relative to the center port) . 

Combination of equations (S) and ( 6 ) yields 


dp 23 2 dp 12 

cos 20 = r + -T 

b 23 q bi 2 q 
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Upon dividing equation (5) with (7) , one obtains 

A P23 


or 


tan 26 = 


/ b 23\ 

A P23 + 2 ' v bl^ / ) A Pl2 


0 = y tan "* 1 


A P23 


a P 23 + 2b r &pi2 


C8) 


where b^ = b 23/bi2* 


UNSTEADY FLOW 


Consider the X component of the Navier-Stokes equations 


3u 3u 3u 3u 

t— + u 7 — ■+ v - 5 — + w 
3t 3x 3y 3z 


3h 1_ 3p_ y. / 3 2 u + 3 2 u . 3 2 u 
3x - p 3x ‘ p ^ 3x2 8y2 * 3z2 


( 9 ) 


For irrotational flow 


3w _ 3v_ 

3y " 3z 1 


3u 

3z 


3w 

3x 


3 v 
3x 


dU 

3y 


hence, equation (9) can be written 


lii 3_ 

( u 2 + V 2 + w 2 \ 

_ 3h 1 3fi. u _3_ 

( 3u_ 3v_ 3w \ 

3t 3x 

\ 2 ) 

. 3x ‘ p 3x p 3x 

y3x + 3y + 3z J 


Since for incompressible f 1 ow 


3u . 3v , 3w . 

vrr T -^77 + — - u 

3x &y op 


one obtains 
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3u + J_ 


at ax 

V 2 , 

where V 2 = u 2 -s 

= V 2 

Similarly 

for 

av _a_ 

( V2 

at + ay 


aw a 

h_L 

at + az 

V 2 


Similarly for the Y and Z components one obtains 


♦*£♦*$*« 


For potential flow 


3 if) 

u " 3X » V “ 3y ' ” ~ Ji 

the time derivatives of the velocity components become 


3u _ 3 2 d> 3v _ 3^4> 3w _ 3 2 d) 

at " 3t3x ’ at atay ’ at ataz 


Therefore equation (10) may be written 


ax l at / ax l 2 s P ) 


(10a) 


(10b) 

(10c) 


_a_ /ajA + _a_ [ 
ay l at J ay \ 2 


a v , p 
+ -rr: . -=r- + gh + 


Cllb) 


az at az \ 2 


J_ li + e h + E B 


(11c) 


Ignoring gravitational effects, one obtains upon integration 

+ ii* E = pm 
at 2 p ' rlZ) 
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where F(t) is a function of time only. Since the above relation holds for 
the entire field of motion at a given time t = tj we may write for ports 2 
and 3 


fit . vL.e\ 

\3t 2 p J 


= f|i , vf. , e) 

\3t 2a . 

at point 2 \ 'at point 


hence 


D vi * t] - [i * 4] ■ ■ (# l - (If) 


Thus, based on the equation of motion the simultaneous pressure difference 
between ports 2 and 3 becomes 


P2 - P3 i ? 2 

- 7 — = J cvi - V 2 ) + 


op 

3t 


or 


Ap2 3 


10 2 

J PCVi - v 2 ) 




(15) 


The first term on the right-hand side of equation (13) is independent of time, 
while the second term may depend on both space and time. It may well be that 
the first term could be replaced by equation (5) , while the second term could 
be replaced by some function of the kinetic energy, space, and time. Thus 


a P2 


3 = — PU 2 b 2 3 sin 29 + j pU 2 • a(x,t) 


(14) 


Since the probe performs an oscillating motion, the space and time variables 
may be interchanged with velocity and frequency of oscillations . Introducing 


e 23(U>f) “ 


a(x,t) 

b 23 sin 29 (t) 
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upon substitution into equation (14) one obtains 


AP2 3 “ 3 ^2 3 


1 + e 23 (U,f) 


sin 20 (t) 


where q = pU 2 . Similarly 


Api 2 = q b 12 — 

/2 


i[" 


ei 2 (U 5 f) 


sin - 20 


(IS) 


(16) 


DETERMINATION OF b 23 , b 12 , and e(U,f) 

l. 

The constants b 23 and b 32 can he determined from experimental data 
obtained from wind tunnel tests. Equations (5) and (6) can be rewritten as 


and 


b 23 = 


A P23 

q sin 26 


(17a) 


Api2 


>12 


t q si„(f - 26) 


(17b) 


Since all the terms on the right-hand side of equation (17) can be measured in 
steady flow tests, the sphere constants b 23 and b 32 can be evaluated 
(see figures 2b and 2c) . 

The unsteady effect on e(U,f) was expected to be strongly dependent on 
the oscillating frequency. The probe pitching angle 0(t) can be written as 

0(t) = 0 Q + a sin(27rft) (18) 

where 0^ is the reference displacement angle, a is half of the amplitude 
angle, and £ is the oscillating frequency. By comparison of equations (5) 


and (6) and (15) and (16), it appears that for the case f = 0, e ia (l),£) = 0, 
and e 23 (U,f) = 0. 

In order to study ej 2 (U,f) and e 23 (U,f) rewrite equations (15) and 
(16) as 


S- * 29 1«) 

and 

ijf-l 1 * • sin (26 It)) 


(19) 


( 20 ) 


The difference between unsteady and steady motion then can be analysed, and 
the value for e 2 3 and e 12 can thus be found. These studies axe still 
being pursued at the time of writing this report. 
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of oscillating probe. 
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Figure 5(b). Velocity vectors during upstroke and downstroke 
at the same 0 angle. 
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(a) Variation of nondimensional probe output Ap/q with pitch 
angle in the range of +30 to -30 degrees 


Figure -t. Probe calibration under steady flow conditions, ODU probe 
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Figure 4. Continued. 
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(d) Variation of ratio b r with Reynolds number 
Figure 4. Concluded 







(a) Ap2 3/ »| at Reynolds no. 3.04 x K)‘* 

I-' i gure (>. Probe characteristics under oscillating conditions lor various frequencies when 
displacement angle 6 ^ = 0.6 degrees. 
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Pitch angle, <• . degrees 
(c) Ap^/q at Reynolds No. 3.0 * lO 1 ' 
Figure 6. Continued 
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(c) A p [ 2 /‘| ‘‘t Reynolds No. 2.38 x 10 
I’igure 7. Continued. 












(I>) AP 23 /M at Reynolds No. (>.25 * H> ' 
ligure 8. Continued. 
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Figure 9. 






(a) Variation of AP 23 /M 

figure 10. Variation of pressure differential with frequency during a period T for constant 
Reynolds number, lie = 0.25 x 1 0 1 * . 
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Normalized time of oscillations 
(l>) Variation of Apj 2 /q 


Figure 10. Continued. 
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